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Abstract. When finding the nonzero eigenvalues for Hamiltonian eigenvalue problems it is especially 
important to locate not only the unstable eigenvalues (i.e., those with positive real part), but those which 
are purely imaginary but have negative Krein signature. These latter eigenvalues have the property that 
they can become unstable upon collision with other purely imaginary eigenvalues, i.e., they are a nec- 
essary building block in the mechanism leading to the so-called Hamiltonian-Hopf bifurcation. In this 
paper we review a general theory for constructing a meromorphic matrix-valued function, the so-called 
Krein matrix, which has the property of not only locating the unstable eigenvalues, but also those with 
negative Krein signature. These eigenvalues are realized as zeros of the determinant. The resulting finite 
dimensional problem obtained by setting the determinant of the Krein matrix to zero presents a valuable 
simplification that can be solved with substantially less computational overhead than the full diagonal- 
ization of the linearization matrix. In this paper the usefulness of the technique is illustrated through 
prototypical examples of spectral analysis of states that have arisen in recent experimental and theoreti- 
cal studies of atomic Bose-Einstein condensates. In particular, we consider one-dimensional settings (the 
cigar trap) possessing real-valued multi-dark-soliton solutions, and two-dimensional settings (the pancake 
trap) admitting complex multi-vortex stationary waveforms. 
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1 . Introduction 

Hamiltonian eigenvalue problems have a time-honored history, as they arise in numerous applications 
stemming from fluid mechanics, celestial mechanics, optical and atomic physics among many other disci- 
plines; see for some recent examples the books [17, 26, 32]. Especially in higher dimensional settings these 
problems can rapidly become fairly computationally intractable, at least as concerns providing the full 
diagonalization of the relevant matrix (e.g., when it stems from the linearization around two-dimensional 
vortex structures or three-dimensional vortex-rings [20]). It is therefore highly desirable to be able to re- 
duce the computational overhead of this procedure, by providing a technique that can capture the main 
features of the linearization spectrum through suitable reductions to a finite dimensional eigenvalue prob- 
lem. It is the aim of the present work to provide a general overview, as well as a systematic set of case 
examples of such a method. The approach that will be developed will be based on the so-called Krein 
matrix [11]. The Krein matrix is a meromorphic matrix -valued function constructed on the basis of a refor- 
mulation and subsequent projection of the relevant eigenvalue problem to a finite dimensional subspace. 
The construction is such that the eigenvalues are realized as points for which the Krein matrix is singular. 
Hence, the computation and visualization of the determinant of this matrix can serve as a tool to identify 
the eigenvalues of the original problem. 

Our presentation will be structured as follows. In Section 2, we first provide an overview of the Hamiltonian- 
Krein index theory for Hamiltonian eigenvalue problems. Afterwards, we focus on the construction of the 
Krein matrix and provide a summary of its properties. In Section 3 we consider specific examples stem- 
ming from the application of the Krein matrix analysis to vortex dynamical states that are of intense re- 
cent interest in the field of atomic physics. In particular, we consider single (unit charge) vortex states 
that are presently fairly routine to produce/observe [18], but also examine the case of the recently stud- 
ied vortex dipoles. The latter were not only produced by dynamical experiments involving flow past an 
obstacle [25], or quenching through the Bose-Einstein condensation transition [5], but their dynamical 
properties and structural robustness were studied at considerable length in recent works in the physical 
literature [21, 24, 28, 30]. For reasons of completeness, we also examine in Section 4 one-dimensional 
analogs of such vortex states; namely, multi-dark soliton structures that have been studied theoretically 
(see e.g. [2, 4] and references therein) and also have been observed experimentally in [29, 31]. Finally, 
in Section 5 we summarize our findings and present our conclusions, as well as some potential topics for 
future study. 

Acknowledgments. TK gratefully acknowledges the support of the Jack and Lois Kuipers Applied Mathe- 
matics Endowment, a Calvin Research Fellowship, and the National Science Foundation under grant DMS- 
1108783. PGK and DY gratefully acknowledge support from US AFOSR under grant FA9550-12-1-0332. 
PGK also acknowledges support from the NSF under grants DMS-0806762 and CMMI-1 000337 and also 
from the Alexander von Humboldt Foundation and the Binational Science Foundation. 



2. Hamiltonian spectral theory (a review) 

Consider the Hamiltonian eigenvalue problem 

J£u = Xu, (2.1) 

acting on a Hilbert space X with inner-product (•, •). The operator J is skew-symmetric bounded operator 
with a bounded inverse, and the operator C : Y i— » X is Hermitian with a compact resolvent. The space Y C X 
is assumed to be dense. Under these assumptions it is well-known that the spectra of JC, namely o{JC), is 
all point spectra, each eigenvalue has finite multiplicity, and infinity is the only possible accumulation point 
of the eigenvalues. If we further assume that each of the operators has zero imaginary part, i.e., Im( J') = 
Im(£) = 0, then the eigenvalues satisfy the quartet symmetry that if A e a{JC), then the set {+A,+A| c 
a(JC). Furthermore, the algebraic multiplicities of each of the eigenvalues in the quartet matches, e.g., 
ma(A) = ma(A). 



T. Kapitula, P. Kevrekidis, and D. Yan 



3 



2.1. The Hamiltonian-Krein index 

The Hamiltonian-Krein index theory has a long history (see [1, 3, 9, 14, 15, 27] and the references therein). 
The index theory is used to related a(C) - in particular, n{C) - to the total number of A e a(J C) with positive 
real part. 

The details for the following discussion can be found in, e.g., [3], and an abbreviated version is in- 
cluded here for the sake of completeness. The eigenvalue problem (2.1) arises from linearizing a particular 
Hamiltonian system about some type of steady state solution (solitary wave, spatially periodic wave, etc.). 
The underlying system has N symmetries, which means that dim[ker(£)] > N, as each symmetry generates 
an element of the kernel, and the kernel elements generated in this fashion are linearly independent. We 
will henceforth assume that dim[ker(£)] = N, with ker(£) = span{(j()i,...,(j()2vl- Now, the generalized kernel 
is found by solving Cii = J~^(f), where (p e ker(£). Now, the symmetries generate conserved quantities, 
and it turns out to be the case that the manner in which these quantities are generated allows us to solve 
£u = J~^<pj for each ; = 1,...,JV: denote these solutions as xfjj. Thus, mg(0) = JV, and ma(0) > 2JV. If we set 
D 6 ]rN><^^ as Dij = (^piXipj), then by the Fredholm alternative it will be the case that ma(0) = 2JV if and 
only if D is nonsingular: this will henceforth be assumed. 

In tlx is a solution to (2.1) for X ^ 0, then by the Fredholm alternative it must be the case that for any 
(p e ker{C), 

= -(J-l„^^,(/,) = {ux,J-^<p) = {ux,C^). 

Consequently, the eigenvalue problem is not solved on all of Y, but is instead solved on the N co-dimensional 
constrained space S-*-, where S = STpan{Cipi, . . . , Cipj^]. Thus, it will be important not to calculate n(£), but 
instead n(£sx), where Ljx = Pg±£Ps± : S-"- i-^ S-"-, and Pjx : X S-"- is the orthogonal projection. It was most 
recently shown in [12] that 

n(£s^)=n(£)-n(D); (2.2) 

thus, the matrix D also plays a significant role in determining the number of negative eigenvalues for the 
constrained operator. 

We are now ready to state the main result. Let denote the total number of positive real-valued 
eigenvalues (including multiplicity), and k^. the total number of complex eigenvalues with positive real and 
imaginary part (again, including multiplicity). Under our assumptions it will be the case that k^. is an even 
integer. For each nonzero eigenvalue A e ilR^ with associated eigenspace Ex, let 

kr{X) :=n(CE,). 

Here we are using the notation that for a self-adjoint operator T and a subspace Z with basis [zi,...,z^], 

The quantity k7{X) is known as the negative Krein index of the eigenvalue, and if k7{X) > 1 the eigenvalue 
is said to have negative Krein signature. Although we will not need the result, it can easily be shown that 
/cj~(A) = A;j~(A). The total negative Krein index is given by 

Aea(J/:)n(iIR+\|0)) 

The Hamiltonian-Krein index is the weighted sum of these indices; namely; 

Ku,m = K + 2k, + 2k[. (2.3) 

The doubling of the indices takes into account the conjugated eigenvalues A which appear in the closed 
fourth quadrant of the complex plane. The major result is that the Hamiltonian-Krein index is intimately 
related to the number of negative eigenvalues of the constrained operator, i.e., 

Kii,^ = n(C)-n{D). (2.4) 

Remark 2.1. A few of implications of (2.4) are as follows: 
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(a) if ^Ham is odd, then k^>l,so that the underlying wave is spectrally unstable 

(b) if -Knam = 0, then it is generically the case that the wave is orbitally stable 

(c) we see that there is a finite and prescribed number of (potentially) unstable eigenvalues, which, 
as we will see in subsequent sections, greatly assists in a numerical search for these eigenvalues. 

Remark 2.2. The Krein signature has important implications beyond what is present in the Hamiltonian- 
Krein index. Namely, if two purely imaginary eigenvalues collide, then after the collision they can attain 
a nonzero real part (this is the so-called Hamiltonian-Hopf bifurcation) if and only if they have opposite 
signature. If the two eigenvalues have the same signature, their collision cannot give rise to instability. 

Before we can construct a function whose zeros correspond to eigenvalues, we must first find a self- 
adjoint pencil which is equivalent to (2.1). First suppose that Re A 0. Upon setting v = J~^u, and setting 

it is not difficult to see that solving (2.1) is equivalent to solving 

L^u = Xv, C_v = -Am. (2.5) 

We now write down an equivalent eigenvalue problem for nonzero A. We have that 

ker(£+) = ker(£), ker(C_) = span|J-Vi.---. J"Vn} = J"^ ker(£) : 

the discussion in the previous subsection allows us to say that ker(£+) ± ker(C_). Upon setting Fl : X 
[ker(>C^) © ker(>C_)]-^ to be the orthogonal projection, it is known from [14, Section 3] that for nonzero 
eigenvalues (2.5) is equivalent to the system 

- n£+nu = Av, n£_nv = Xu. (2.6) 

Each of the operators n>C+n are self-adjoint, and the assumption that D is nonsingular implies that each is 
also nonsingular on the range of FT. This allows us to introduce the invertible self-adjoint operators 

n := nc+n, s-^ -.= nc_n, 

and rewrite (2.6) as the quadratic pencil 

{n + x^s)u = 0. 

Similarly, if we now assume that A e ilR""", then the equivalent pencil is 

{TZ + X^)[lmu] = 0, Rew = {lmX)JS[lmu]. 

In conclusion, we have that solving the original eigenvalue problem (2.1) is equivalent to solving the 
linear pencil 

{TZ-zS)u = 0, z:=-X^ (-7r/2 < arg A < Tr/2). (2.7) 

The effect of the eigenvalue mapping is illustrated in Figure 1. Eigenvalues with positive real part and 
nonzero imaginary part are mapped in a one-to-one fashion to eigenvalues with nonzero imaginary part, 
and the four-fold symmetry is reduced to the reflection symmetry {z,zj. The system (2.6) has an unstable 
eigenvalue A with positive real part if and only if the system (2.7) has an eigenvalue z with z < or with 
Imz * 0. 

Let us conclude by computing the Hamiltonian-Krein index for the pencil (2.7). First consider the 
purely imaginary eigenvalues. It is straightforward to show that 

J^E, = 27^£^, 

where z = -A^ e E^, and = ImE^. Since = ReE;^ ±ilm£^, it is consequently the case that 

kr(z) = k.{X) + kr(-X) = 2kr(X), 
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where for z e IR"*" the Krein index is found by computing 

fc^(z) = n(7^£J. 

Thus, the total negative Krein index is the same for the pencil as it is for the original problem. Now consider 
those eigenvalues for the original problem with nonzero real part. Set k^iz) to be the multiplicity of the real- 
valued eigenvalue z e IR""", and let k^{z) be the multiplicity of the eigenvalue for Imz ^ (it is clearly the case 
that /cc(z) = kf-{z)). Since 

{TZ-{±XfS)u^x = 0, 

we clearly have that 

k,{z) = 2k,{X), fcjz) = 2fc,(A). 

Here we are using the notation that k^{X) is the multiplicity of the positive real-valued eigenvalue for (2.1), 
and k^(A) is the multiplicity of an eigenvalue with positive real and imaginary parts. Upon summing over 
all of the eigenvalues, and using (2.3), we see that the Hamiltonian-Krein index for the original problem is 
related to that of the pencil (2.7) in the following manner: 

Lemma 2.3. Consider the linear pencil (2.7) as derived from the eigenvalue problem (2.1). For the pencil let 
k^ denote the total number of negative real-valued eigenvalues (counting multiplicity) , k^- the total number of 
eigenvalues with positive imaginary part (counting multiplicity), and k^ the total negative Krein index of all the 
positive real-valued eigenvalues, where the index for a single eigenvalue z e IR"*" with associated eigenspace is 
given by 

/c^(z) = n(7^£J. 

Then with the Hamiltonian-Krein index given as in (2.4), the eigenvalues for the pencil satisfy 

k, + 2fc, -I- 2k[ = 2J<:Ham- 

We conclude by relating the Hamiltonian-Krein index to n(iS~^) = n(7^) (the equality follows from the 
fact that J has bounded inverse). As for the number of negative directions for S~^, we have 

n(5-l) = n(-( J£J)ker(£)x) = n(/;[j-lker(£)]x) = n(£) -n(D) = ^Ham- (2.8) 

The first equality follows from the definition of and the fact that constrained operator maps the sub- 
space ker(£)^ to itself, the second equality follows from a simple change of variables, the third equality 
follows from (2.2), and the fourth equality is from (2.4). Thus, since the number of negative directions is 
invariant under inversion, with respect to the pencil alone we can rewrite the conclusion of Lemma 2.3 as 

k, + 2k^ + 2kr = 2n{S). (2.9) 

It can be concluded that if the underlying wave is orbitally stable, then S is positive definite; otherwise, the 
operator must be indefinite, although it will necessarily have only a finite number of negative directions. 

2.2. The Krein matrix 

We now turn to the problem of constructing a meromorphic matrix-valued function, the Krein matrix, 
which has the property that it is singular precisely for those values of z which correspond to nonzero 
eigenvalues for the pencil (2.7). The Krein matrix was introduced in [11], and the interested reader should 
consult that work for the details associated with the following discussion (also see [12, Section 3]). 

We now proceed to construct a meromorphic matrix-valued function, the Krein matrix, which will 
be used to find the eigenvalues for the pencil (2.7). We construct the Krein matrix by projecting off the 
negative space of the operator S, which as we have already seen has dimension JCnam- Let N(S) denote the 
^^Ham "dimensional negative subspace of S, and let P : X i-^ N(5)-^ be the orthogonal projection. Define the 
constrained operators 

7?.2 := PT^P, S+ := PSP. 
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Figure 1: (color online) Six sets of eigenvalues and their images under the map. The (red) circles 
denote two quads of complex eigenvalues under the four-fold symmetry j+A, ±A) and their images, 
k(- = 2. The (green) crosses denote two pairs of real eigenvalues j+Aj and their images on the negative 
real axis, = 2. The (blue) boxes denote two pairs of purely imaginary eigenvalues j+A) and their 
images on the positive real axis. The filled square has a positive Krein signature, while the empty 
square has a negative Krein signature, so that fcj" = 1 . 

The operator is positive definite and symmetric, whereas from the Index Theorem [12, Theorem 2.1] the 
symmetric operator 7?-2 satisfies 

n(TZ2) = n(7e) -n(7e-i5,) = ^CRam " nlTe'i^,). 
Finally, since 5+ is positive definite we can define the conjugated operator 

n := s;^^^TZ2S;^^^, 

and note that it satisfies 

n(^) = n(7^2). 

Upon w^riting a potential eigenfunction to the linear pencil as 

11= CjSj + s+, s+ e N{S)-^, 

where {si,...,Sf^^_^^} is a basis for N(S), the eigenvalue problem is equivalent to 

K{z)c = 0, JC(z):=71n(5)-z5n,5)-(^-z)-Jv,p^^^^j. (2.10) 
Here we are using the notation 

S;^^^PTZN{S) := {S^'^^^PTZs- : s" e N{S)}. 

Now, by construction z is an eigenvalue for the pencil (2.7) with corresponding eigenfunction u if and only 
if K{z)c = 0, where c = (ci,...,Cxjj^^)^. If for a particular eigenvalue z it is true that ii & N(S)-'-, then it 
is necessarily true that detK'(z) = 0. On the other hand, if m e N(S)-'-, then c = 0, so it can be the case 
that detK'(z) ^ 0. Now, if Imz ^ 0, then it will always be the case that z is an eigenvalue if and only if 
detiiC(z) = 0. On the other hand, if z e IR~ is an eigenvalue with detiiC(z) ^ 0, then it is the case that for the 
Krein matrix constructed by projecting off of the negative directions of TZ, say Kii(z), we would necessarily 
have det Kn(z) = 0. Finally, if z e IR"*" is an eigenvalue with det K'(z) ^ 0, then the eigenvalue has positive 
Krein signature. In this case the eigenvalue is realized as a removable singularity of the Krein matrix, i.e.. 
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z = Zp is a pole of the Krein matrix for which the residue is the zero matrix. If the eigenvalue has negative 
Krein signature, it will necessarily be the case that detiiC(z) = 0. 

In conclusion, we can use detiiC(z) as a meromorphic function whose zeros correspond to eigenvalues. 
The (potential) singularities of the Krein matrix arise through (TZ - z)~^ ,„ at the eigenvalues of the 

self-adjoint operator TZ. In order to use the Krein matrix to say something about the Krein signature of a real 
positive eigenvalue, it will be helpful to factor the determinant as a finite product. One can easily observe 
that K(z) is symmetric for all z, i.e., K(z)^ = K{z); in particular, it is Hermitian for z e IR. This allows us for 
z 6 IR to extract the Ky^^^ eigenvalues of the Krein matrix, ry(z), hereafter called the Krein eigenvalues, in 
a meromorphic fashion. Thus, instead of finding eigenvalues by looking for the zeros of the determinant 
of the Krein matrix, we can look for the zeros of each individual Krein eigenvalue. There will be precisely 
^Ham of these Krein eigenvalues. 



Figure 2: (color online) This is a cartoon of the graphs of the Krein eigenvalues in the case that 
KHam = 2. The eigenvalues with positive Krein index are denoted by (blue) squares, and those with 
negative index are shown as (red) circles. The dashed vertical lines correspond to vertical asymptotes 
for each Krein eigenvalue. The arrows connect adjacent zeros of Krein eigenvalues. 

For real-valued z, the properties of the Krein eigenvalues are as follows. First, they can be enumerated 
so that 

lim ^>0, 

z— >-oo 2 

SO that each Krein eigenvalue is negative for large negative z. This follows from the fact that <Sjv(5) is a 
negative definite matrix. Second, if Zp is a pole of the Krein eigenvalue rj(z), it will then be the case that 

lim rj(z) = +00. 

Furthermore, if z* is a simple eigenvalue of TZ, so that z* is a pole of the Krein matrix, then it will be the case 
that z* is a simple pole for only one of the Krein eigenvalues. In other words, for all but one of the Krein 
eigenvalues the pole z* of the Krein matrix is a removable singularity. Finally, if Zq is a simple positive zero 
of the Krein eigenvalue rj(z), then it is true that 

Ki^o) = -sign[r^'(zo)]; 

in other words, the slope of the Krein eigenvalue at a zero gives definitive information regarding the Krein 
index of the eigenvalue. This discussion is summarized in Figure 2. 



3. Application: spectral analysis for vortices of the GP equation 



We now wish to use the Krein matrix to identify the eigenvalues of (2.1) which have nonzero real part, or 
which are purely imaginary and have negative Krein signature. The Hamiltonian-Krein index of (2.9) tells 
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us how many of these zeros for the pencil we need to find in order to fully capture all of the (potentially) 
unstable eigenvalues. The real eigenvalues for the linearized problem correspond to negative real eigen- 
values for the pencil, and the eigenvalues with nonzero real part correspond to eigenvalues with nonzero 
imaginary part for the pencil. Regarding the eigenvalues with negative signature, from the discussion of 
the previous section this means that we need to identify the real and positive eigenvalues for the pencil 
(2.7) for which a Krein eigenvalue has a zero, and the slope of the curve at the zero is positive. 

The model under consideration for the case of pancake-shaped Bose-Einstein condensates [5, 18, 24] is 
the (2-i-l)-dimensional Gross-Pitaevskii equation in the dimensionless form 

1 , 
idfU = -—Am + V(x,y)u + \u\ u- fiu. (3.1) 

Here, u is the macroscopic wave function, V{x,y) = Q^(x^ -i-y^j/2 is the external harmonic potential, Q is 
the frequency of the trap, ^ is the chemical potential, and A is the Laplace operator, i.e., A = dxx + dyy For 
the problem discussed herein, we assume Q = 0.2. We begin by assuming that the steady-state problem 
is solved, and we will denote that solution by UQ(x,y) = UQ{x,y) + \ VQ(x,y). Here Uq,Vq are real-valued 
functions. Abusing notation a bit and writing 

u Uq + Wq + e{u +w), < e «; 1, 

we see that at 0{e) we have the linear problem 

dtU = JCu, u = {u,vj^, 

where 

j^, 1\ -^A-ii+V{r) + 3U^ + V^ 2UoVo 



1 /' " \ IUqVq -^A-iA + Vir) + U^ + 3V^ 

Here we take advantage of the fact that the trap is radially symmetric by writing = x^+y^. The eigenvalue 
problem of the form (2.1), i.e., 

JCu = Am, (3.2) 

arises upon using the separation of variables ansatz u i-^ Me'^^ 

With respect to the standard inner-product on L^(]R^) x L^(]R^) the operator J is boundedly invertible 
and skew-symmetric, while the operator £ is self-adjoint. Furthermore, since the potential V{r) grows 
quadratically, and the solution Mo(x,y) decays sufficiently fast as r — > +oo, it is the case that £ has a compact 
resolvent. Following (2.5), we set 

= £, £_ = -JCJ. 

In order to construct the Krein matrix, we must consider the spectral problem on the space perpendicular 
to the kernel of both operators. The fact that solutions to (3.1) are invariant under u(x,y) ii{x,y)e"l' and 
the spatial rotation ii(x,y) u(xcos6 -ysmd,xsm6 +ycos6) means that (generically) the kernel of £ will 
be two-dimensional. 



(3.3) 



ker(£+) = ker(£) = span{(-Vo; Uq)^, (dgpcos(p - ^^(/jVo'O'epsinc/) -l- dgcjyUQ)^} 
ker(£_) = ker(- J£ J) = span{([/o, Vq)^, {d0psm(p + dQ(j)Uo,-dgpcos(j) + decpVof]. 

For the second eigenfunction we are writing the underlying wave as 

Uo{x,y) + iVo(x,y) = p(r,0)e'*<''^', tan0 = y/x. 

If the solution has a radially symmetric density, i.e., p = p(r), then the dimension of each kernel is one with 

ker(£^) = span|(-Vo, Uof], ker(£_) = span{([/o, Vof]. 

With this information at hand we can now compute the restricted operators TZ,S for the pencil (2.7). 
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Figure 3: The single vortex state for a trap strength Q = 0.2 and chemical potential /< = 3. The 
left panel shows the contour plot of the density, and the right panel is the phase plot for the wave 
function. In this case, the total number of grid points isn = 120, and the spatial step size is Ax = 0.25. 

3.1. Single vortex state 



Here we assume that the solution is a vortex of charge one. In other words, the density profile p = y + 
is radially symmetric v\rith p(0) = 0, while the phase rotates by 2n around the vortex core, justifying the 
topological charge of the structure [18, 19]. The density and phase profiles for the single vortex state are 
shown in Figure 3, and the corresponding linearization spectrum is shown in Figure 4. The dependence 
of the relevant eigenvalues as a function of the canonical parameter of the system, namely the chemical 
potential associated with the number of atoms in the condensate, has been quantified previously; see 
e.g. [19, 23]. We will confirm these findings via the Krein matrix, and showcase particular examples for 
a few representative values of fi. As we saw above, the kernels of the operators £+ each have dimension 
one. For small amplitude vortices where the nonlinear interactions are (almost) negligible it can be shown 
via an analysis similar to that presented in [16] that the Hamiltonian-Krein index satisfies K^ia^^ = 2. By 
Lemma 2.3 we then know that for the pencil 

k, + 2k^ + 2kr = A; (3.4) 

thus, if there are two positive real zeros of the Krein eigenvalues which correspond to an eigenvalue with 
negative Krein index, the rest of the spectrum for the pencil must be positive and purely real. In other 
words, if for the original eigenvalue problem there is a purely imaginary eigenvalue with negative Krein 
index, the rest of the spectrum must be purely imaginary with positive Krein index. Consequently, once we 
numerically find one purely imaginary eigenvalue with negative Krein signature, or one set of eigenvalues 
with nonzero real part which satisfy the Hamiltonian eigenvalue symmetry, we need not search for any 
additional unstable eigenvalues. They simply do not exist. 

Regarding the numerical solution and analysis of the problem, we first discretize the differential opera- 
tor via a centered finite difference scheme using n data points and a spatial interval of Ax. The single vortex 
state - indeed, any of our stationary solutions - is found by applying Newton's method to the discretized 
problem. As for the spectral problem, the smallest (in norm) eigenvalues are computed using the MATLAB 
function "eigs" on the discretized version of the operator 

We now turn to a direct comparison of the numerical results for eigenvalues obtained from the linear 
stability analysis and from the Krein matrix. Here, we have confirmed the above findings with fi = 0.45 in 
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Figure 4: The real and imaginary parts of the eigenfrequencies defined hy co = 
state when n = 120 and Ax = 0.25 as a function of the chemical potential ja. 
stable for all considered values of the chemical potential. 



iX for the single vortex 
The wave is spectrally 



Single vortex n=0.45, Q.=0.2 




0.05 0.1 0.15 
z 



0.01 



Re(X) 



Figure 5: The numerically generated spectral plot for n = 48, fi = 0.45, Ax = 0.45. In the left panel 
the numerically generated plot of the XHam = 2 Krein eigenvalues is shown. For this problem the 
Krein matrix is a meromorphic multiple of the identity; hence, the two Krein eigenvalues coincide. 
The (red) crosses are the poles of the Krein matrix. If the pole is removable (e.g., for z ~ 0.155), then 
it corresponds to an eigenvalue with positive Krein signature. Here we see two positive real zeros 
of the Krein eigenvalues for which the functions have positive slope. Thus, for the pencil = 2 
with fc]- = fcc = 0, which means that for the original problem (3.2) = 1 with k^ = k^ = 0. The wave 
is spectrally stable, but is not a (local) minimizer. In the right figure the eigenvalues for JL are 
denoted by (blue) circles. The (green) crosses represent the poles of the Krein matrix, and the (red) 
triangles are the eigenvalues of JC which are realized as zeros of the Krein eigenvalues. The purely 
imaginary eigenvalue with negative Krein signature is shown as a filled (red) triangle. 
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Figure 6: The numerically generated spectral plot for n = 48, fi = 0.75, dx = 0.45. In the left panel the 
left plot is the numerically generated plot of the JCnam = 2 Krein eigenvalues. For this problem the 
Krein matrix is a meromorphic multiple of the identity; hence, the two Krein eigenvalues coincide. 
The (red) crosses are the poles of the Krein matrix, and the removable singularities (e.g., z ~ 0.08) 
correspond to eigenvalues with positive Krein signature. Here we again see two positive real zeros of 
the Krein eigenvalues for which the functions have positive slope. The wave is then spectrally stable 
with k~ = 1 for the problem (3.2). The notation used in the right figure is similar to that in Figure 5. 



Figure 5, and = 0.75 in Figure 6. In both of these figures a spatial discretization spacing of Ax = 0.5 was 
used. Similar results were found with smaller values of Ax. In the right panel of each figure the spectrum of 
J'C was computed for both the original formulation of (3.2) and the corresponding pencil formulation (2.7), 
and the results of the two were found to agree within the accuracy of the eigenvalue solver. These spectral 
results were then compared with the prediction of the Krein matrix in the right panel. The agreement is 
excellent. 

Unexpectedly, it turns out to be the case that for the problem at hand the Krein matrix is a meromorphic 
multiple of the identity; hence, the two Krein eigenvalues coincide. In both figures it is clear that there are 
then two zeros of the Krein eigenvalues at which the slope of the curve is positive: this corresponds to a 
purely imaginary eigenvalue with negative Krein signature. As stipulated by Hamiltonian-Krein index in 
(3.4), the rest of the spectrum must be purely imaginary, and this is indeed seen to be the case. 

In terms of /^-independent eigenvalues within the spectrum, there is not only the eigenvalue at A = 
which is due to the phase U(l) gauge invariance of the model, but also a double eigenvalue at A = +0.2i. 
The latter frequency of double multiplicity is the so-called dipolar mode of the condensate and pertains to 
a symmetry (oscillation of the entire condensate cloud in the x- or y- direction with the trap frequency), 
and is hence invariant with respect to variations in j^. One of these modes pertains to a pole of the Krein 
matrix, while the other is a zero of a Krein eigenvalue. The eigenvalue with negative Krein signature always 
lies between the origin and this double pair and is known to tend to the origin as the chemical potential fi 
increases [19, 23]. Since the eigenvalue has negative signature, as predicted by the theory it is realized as a 
zero of a Krein eigenvalue. 

To showcase the potential relevance of poles, the residues of the poles of the Krein matrix are computed. 
If the numerically computed residue of the pole is of 0(10"^^), then we say that the pole is removable, and 
hence it corresponds to an eigenvalue with positive Krein signature. The residue is computed via the 
numerical integration 

If 1 " 

Res(K(z),Zp) = — ^^K{z)dz ~ — ^li:(zp,,)Azp,,. 

Here Zp is the relevant pole of the Krein matrix, ti is the number of integration points on the simple pos- 
itively oriented closed contour C which surrounds the pole (here the contour is without loss of generality 
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Figure 7: The vortex dipole state for a trap strength Q = 0.2 and chemical potential ja = 3. The 
left panel shows the contour plot of the density, and the right panel is the phase plot for the wave 
function. In this case, the total number of grid points is n = 120, and the spatial step size is dx = 0.25. 



chosen to be a square centered on the pole), Zp is the point on the integration contour around the pole, 
and Azp , is the segment on the integration path. 

When the chemical potential = 0.45 the first relevant pole we choose is located at Zp = 0.0021. For the 



increment of |Azp , | = 5*10 ^ it is seen that 

-8.4* 10-12 -12.1 * IQ-^^ -4.5* 10-22 -i9.5* 10-20 

ic 



i 



^^^^^^-[ -2;8* 10-21-11.1*10-1'^ -8.4*10-12-12.1*10-14 



so that I Res(iiC(z), 0.0021)1 = 0(10-12). The simple pole is removable, and corresponds to a purely imaginary 
eigenvalue (the first red cross above zero in Figure 5). The fact that the singularity is removable is evidenced 
by the fact that neither of the Krein eigenvalues has a singularity at the point. Now consider the pole at 
Zp = 0.0478. It is seen that 

r „, . , _/ -10.0057 -8.6*10-1^-13.0*10-1^ 

y^/((z)dz-^ _55,io-i«-i2.2*10-i5 -10.0057 

so that I Res(iiC(z), 0.0478)1 = 0(10-^), which is nonzero by our criterion. We thus see that the pole is not 
removable because the Krein eigenvalues have a singularity at that point. Since the singularity is not 
removable, this point does not correspond to an eigenvalue. 



3.2. Vortex dipole state 

We now turn to a spectral analysis for the vortex dipole state, which is a stationary vortex-antivortex state 
(see Figure 7 for a typical example of its density and phase). When Q = 0.2 the vortex dipole state exists 
for n > 0.68. It is interesting to note here that as shown in [21, 23], such states do not exist at the linear 
limit, but only bifurcate through a supercritical pitchfork (symmetry-breaking) event at a critical point 
from the dark soliton, which corresponds to the real-valued (first) excited state of the two-dimensional 
harmonic oscillator. Since the density is not radially symmetric, it will be the case that dim[ker(£+)] = 2. 
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Figure 8: The real and imaginary parts of the eigenfrequencies defined by w = iA for the vortex 
dipole state when n = 480 and dx = 0.0625 as a function of the chemical potential j^. The wave is 
spectrally unstable with k^- = 1 when Im a) > 0. 



Furthermore, it is seen numerically that n(5) = K^^^ = 2, so that by Lemma 2.3 we have the same index 
count for the pencil as (3.4); namely, 

k, + 2k^ + 2kr = 4. 

As in the previous example, there will be two Krein eigenvalues to be plotted. Unlike the previous 
example, however, the Krein matrix will not be a meromorphic multiple of the identity; hence, the Krein 
eigenvalues will not generically overlap. If there are two positive real zeros of the Krein eigenvalues for 
which the curves have positive slope, then the spectrum for the pencil will be positive and purely real with 
k^ = 2. Otherwise the underlying wave will be unstable. As we will see, in this example this instability 
will arise when k^ = I. It is important to note here that for the pencil it must be the case that if k^ = 2, 
then each Krein eigenvalue has a positive zero at precisely the same point. This fact is a consequence of the 
relationship of the spectrum between the original eigenvalue problem and the pencil; in particular, the fact 
that the spectrum of the pencil "doubles up" the nonzero spectrum for the original problem. 

The spectrum once again features the twofold degenerate dipolar modes which are associated with the 
oscillation with the trap frequency of the whole condensate in the x- and y- directions. In addition to these 
modes, there exists a negative Krein signature mode, which in this case collides with the mode departing 
from zero. This, in turn, leads to the formation of an interval of ;i-values where the spectrum possesses 
complex eigenfrequencies associated with oscillatory instability. Therefore, when studying the computa- 
tion of the Krein matrix of the vortex dipole state and comparing its eigenvalues against the linearization 
analysis, we select two cases, namely ^ = 0.71 and = 0.80. As we see in Figure 8, the former is before the 
collision of the two modes (i.e., k^ = 0), and the latter is after the collision has occurred (i.e., k^ = I). For 
larger values of ^ the complex-valued eigenvalues return to the imaginary axis. 

When fi = 0.71 we show the computation of the Krein eigenvalues (left panel) and the linearization 
spectrum (right panel) in Figure 9. It can be observed that the Krein eigenvalues are very consistent with 
those via linear stability analysis except for the real eigenvalues, which are O(10~^) and are the numerical 
approximation of the known zero eigenvalues. For the original eigenvalue problem the Hamiltonian-Krein 
index is 

K + 2kr + 2kT = 2. 



Since for the pencil it is the case that each of the Krein eigenvalues has a positive zero with positive slope, 
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Figure 9: The numerically generated spectral plot for the vortex dipole when Q = 0.2 and f.i = 0.71. 
The left panel is the plot of the two Krein eigenvalues: the (red) crosses represent the poles of the 
Krein matrix, and the removable singularities correspond to eigenvalues with positive Krein signa- 
ture. Note that fcr = 2 for the pencil, and that the zeros of the Krein eigenvalues which correspond 
to eigenvalues with negative Krein signature coincide. In the right panel the (blue) circles are the 
eigenvalues for JC, the (red) triangles are the eigenvalues which correspond to zeros of the Krein 
eigenvalues, and the (green) crosses are the eigenvalues which correspond to removable singularities 
of the Krein matrix - these eigenvalues have positive Krein signature. The purely imaginary eigen- 
value with negative Krein signature is shown as a filled (red) triangle. The labeling of eigenvalues in 
the right panel is the same as for Figure 5. 
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Figure 10: The numerically generated spectral plot for the vortex dipole when D = 0.2 and = 
0.80. The left panel is the plot of the two Krein eigenvalues: the (red) crosses represent the poles 
of the Krein matrix, and the removable singularities correspond to eigenvalues with positive Krein 
signature. Since there are no positive real zeros with positive slope, it must be the case that for the 
pencil fcc = 2. In the right panel the (blue) circles are the eigenvalues for J L, the (red) triangles are 
the eigenvalues which correspond to zeros of the Krein eigenvalues, and the (green) crosses are the 
eigenvalues which correspond to removable singularities of the Krein matrix - these eigenvalues have 
positive Krein signature. The labeling of eigenvalues in the right panel is the same as for Figure 5. 
The only difference is that the (blue) stars represent the complex eigenvalues with nonzero real part 
as found by the zeros of the Krein eigenvalues. 
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Figure 11: Phase plot for the Krein eigenvalues when = 0.80. The colors represent the argument of 
the Krein eigenvalues, arg(A(£i (z))) and arg(A(£2(z))). The two white spots on each of the two panels 
represent the complex eigenvalues, corresponding to z = 0.021 ±i0.09. 



we have that for the original problem fcr = 1. The eigenvalue with negative Krein signature is denoted 
by a (red) filled triangle in the right panel. As a consequence of the index theory we know that all other 
eigenvalues must be purely imaginary with positive Krein signature. 

When 14 = 0.80 we show the computation of the Krein eigenvalues (left panel) and the linearization 
spectrum (right panel) in Figure 9. It can be observed that the Krein eigenvalues are very consistent with 
those via linear stability analysis except for the real eigenvalues, which are 0(10"^) and are the numerical 
approximation of the known zero eigenvalues. The Krein eigenvalues have no zeros with positive slope; 
hence, it must be the case that for the pencil = 2; numerically, the complex -valued eigenvalue for the 
pencil is = 0.03 + i0.148. This in turn implies that for the original linearization k^- = 1; numerically, the 
eigenvalue is A = +^l-zt. As a consequence of the index theory we know that all other eigenvalues must be 
purely imaginary with positive Krein signature. In order to better illustrate the complex quartet of eigen- 
values, when }i = 0.80 we show the phase plot of each of the Krein eigenvalues in Figure 11. The axes denote 
the complex z-plane, and the colorbar corresponds to the phase of the eigenvalues of the Krein matrix. The 
points where the phase becomes singular correspond to the location of the complex eigenvalues, and are 
thus the points of interest to us within this graph. These points are labeled by the white spots. 

4. Application: spectral analysis for multi-solitons of the GP 

equation 

The model under consideration for this application is the (l-i-l)-dimensional Gross-Pitaevskii equation in 
the dimensionless form 

1 2 
idfU = - — dxxU + V{x)u + \u\ u-fiii, (4.1) 

where now V{x) = V)}x^ll is the external harmonic potential. When doing numerical computations for this 
problem we assume Q = 1. For this problem we are interested in studying the spectrum associated with 
real-valued (multi-soliton) solutions. These will be denoted by Ui^ix), and they will be realized as solutions 
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Figure 12: Linear stability analysis for the 3-dark-soliton states: left panel is Im(A) versus }i, while 
the right panel is the Re(A) versus }i. The first critical value of /i where one pair of the complex 
eigenvalues becomes purely imaginary is approximately n = 4.46. The second critical value where 
the other pair of complex eigenvalues becomes purely imaginary is approximately }i = 10.29. 

of the ODE 

-^^.«Uo + V{x)Uo + Ul - /'Uo = 0. 

The eigenvalue problem associated with the real-valued solution will be exactly that as given in (3.2), 
except that now £ will be diagonal, i.e., C = diag(>C+,£_), with 

A = -\dl - F + ^(^) + 3t/2, £_ = -Idl -n+ V{x) + Ul 

The fact that C is diagonal, i.e., the eigenvalue problem is in canonical form, means that not only can 
more be said about the Hamiltonian-Krein index for the original problem, but the index for the associated 
pencil will also change. This latter amendment follows from the fact that the pencil can be formed directly 
without first passing to the intermediate stage of "doubling-up" the eigenvalue problem. 

The operators £+ are both Sturmian: that, combined with the fact that Uq will decay exponentially fast 
as \x\ +00 implies that the spectrum of each operator will be purely real, and will be composed solely of 
simple eigenvalues. The operator C+ will be invertible, while C-Uq = implies that dim[ker(£_)] = 1. If we 
let n : X spanjUol""" denote the orthogonal (spectral) projection, and set 

7^ = ^£+^, S'^ =nc_n, 

then the search for nonzero eigenvalues for the original problem is equivalent to finding the spectrum of 
the pencil 

{TZ-zS)u = 0, z = -X^ 

(compare with (2.7)). 

The Hamiltonian-Krein index satisfies 

^^Ham=n(7^) + n(5) 

(compare with (2.4)). Since FT is a spectral projection, it will be the case that n{S) = n(£_). As for n{TZ), it is 
the case that 

moreover, upon using the fact that 

C4d^Uo) = Uo ^ {Uo,C-^'Uo)=^d^{Uo,Uo) 
we have the rewritten expression 

n(Tl)=n{C^)-n(d^{Uo,Uo)). 
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Figure 13: The left panel shows the power P(;<) (see (4.2)) versus the chemical potential j^: the 
power plotted here is scaled by the max{P(fi)). The right panel demonstrates the stationary solutions 
as fi varies. The solid blue, dashed green, and dash-dotted red curves represent the cases of |( = 
4.2,4.6,10.3 respectively. Here the grid size is « = 800 and the spatial interval is Ax = 0.025. 



In other words, the slope of the power curve P(}i) = (Uq, Uq) has a direct influence on the number of negative 
directions of the operator C^. In conclusion, the Hamiltonian-Krein index for the canonical problem is 

^^Ham=n(£J + n(/:_)-n(PV))' PW = (UoMo). (4.2) 

There is also an instability criterion. A classical result of [6, 10] (recently reproven in [12] using the Krein 
matrix) allows us to say that a lower bound on k,. is caused by a difference in the number of negative 
directions of TZ and S. In other words, 

K>\n{C_)-[n(C^)-n{P'{^l))]\, (4.3) 

which in particular implies that /c^ > 1 if |n(£_) -n(£^)| > 2. 

Regarding the Hamiltonian-Krein index for the pencil, the fact that the original system no longer needs 
to be "doubled-up" in order to be put into canonical form means that we just need to count eigenvalues 
with respect to the eigenvalue mapping z = -A^) (see Figure 1). In particular, each of the individual indices 
doubles, which in the end means that the index for the pencil is precisely that for the original problem, and 
is given in (4.2). Additionally, it should be noted that in the construction of the Krein matrix the size no 
longer directly depends on -Knam; instead, it will be of size n(£_) = n(S). 

In the previous section we considered the spectral stability of solutions for which two Krein eigenvalues 
could be used to locate the spectrum. We now consider an example for which three Krein eigenvalues are 
needed. The steady-state solution Uq to be considered, hereafter called a 3-dark-soliton, is realized as a 
continuation of the Gauss-Hermite function (8x^ - 12x)e~^^^^ from the chemical potential ja* = 3.5 (recall 
that we assume Q = 1). Let us first determine the Hamiltonian-Krein index associated with this solution. 
A weakly nonlinear analysis along the lines of that presented in, e.g., [13], shows that for fi*\ small, 
the power satisfies the relationship P'(fi) > (the details are left for the interested reader). As we see from 
the numerics (see Figure 13), this relationship holds for all values of fi under consideration. Regarding the 
indices of the operators C+, the combination of Uq having three zeros and Sturm-Liouville theory implies 
that n(£_) = 3; consequently, in the Krein matrix analysis there will be three Krein eigenvalues to follow. 
Regarding the operator C+, it is the case that in the weakly nonlinear limit n(£+) = 3: again, this relationship 
holds for all under consideration. Appealing to (4.2) we see that K^ia^ = 6; unfortunately, the lower bound 
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Figure 14: The spectral picture when ft = 4.2, for which k^- = 2 and fcr = 1. In the left panel is the plot 
of the three Krein eigenvalues. As predicted by the theory, there is only one positive zero of a Krein 
eigenvalue for which the slope is positive. In the right panel the eigenvalues for JC are denoted by 
(blue) circles. The (red) triangles are the eigenvalues of JC which are realized as zeros of the Krein 
eigenvalues. The (blue) stars represent the complex eigenvalues with nonzero real part as found by 
the zeros of the Krein eigenvalues. The eigenvalue with negative Krein signature is marked with a 
(red) filled triangle. Note that this eigenvalue is close to an eigenvalue with positive Krein signature, 
which implies that the system is close to a Hamiltonian-Hopf bifurcation point. 
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Figure 15: Similar to Figure 14, except now fi = 4.46 and = 1 with k. = 2. 



of (4.3) on leads to no definitive conclusion. Regarding the types of (potentially) unstable eigenvalues, 
we have the following possibilities: 
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In Figure 12 we show the location of the first few eigenvalues from the linear stability analysis as a 
function of }i. It is always the case that = 0, while /c^ decreases from two to one to zero. The decreasing 
of fcc implies that fcr increases from one to two to three. The critical values of for which fc^ decreases 
are approximately }i ~ 4.465 and }i ~ 10.29. In Figure 14 we show the plot of the Krein eigenvalues (left 
panel) and full linearization spectrum (right panel) when p. = 4.2, when k^. = 2 and = 1. The eigenvalue 
in the upper half-plane with negative Krein signature is denoted with a (red) cross. In Figure 15 we show 
the plot of the Krein eigenvalues (left panel) and full linearization spectrum (right panel) when ^ = 4.46, 
when k^ = I and = 2. The two eigenvalues in the upper half-plane with negative Krein signature are 
denoted with (red) crosses. Finally, in Figure 16 we show the plot of the Krein eigenvalues (left panel) and 
full linearization spectrum (right panel) when = 10.30, when k^ = 3. The three eigenvalues in the upper 
half-plane with negative Krein signature are denoted with (red) crosses. 

Remark 4.1. In general, when considering N-dark-solitons to (4.1) it will be the case that P'(fi) > with 
n(C^) = n(C_) = N. Consequently, when looking for the 2N (potentially) unstable eigenvalues for the 
linearized problem it will be the case that N Krein eigenvalues must be plotted. 



5. Conclusion 

In the present work we have revisited Hamiltonian skew-adjoint eigenvalue problems that typically arise 
in the linearization around a stationary state of a Hamiltonian nonlinear partial differential equation. We 
presented a brief overview of the known facts for the eigenvalue counts of the corresponding (potentially) 
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Figure 16: Similar to Figure 14, except now ja = 10.30 and fcr = 3. 

unstable spectra. We especially focused on a novel, but straightforward, plan to implement finite dimen- 
sional techniques to locate this spectrum via the singular points of the meromorphic Krein matrix. We 
illustrate the value of the approach by considering realistic problems for recently observed experimen- 
tally multi-vortex and multi-soliton solutions in atomic Bose-Einstein condensates. We believe that this 
approach can provide a valuable alternative to the highly-demanding computations that require a diago- 
nalization of the linearization matrix, especially in two- and three-dimensional settings. 

Naturally, there are numerous possibilities for further exploration. It would be especially relevant for 
the atomic physics applications to use the present method to examine the details of the spectra of three- 
dimensional solutions such as vortex rings [20]. On the other hand, from a methodological perspective, 
it would be also relevant to consider this approach for the case of solutions of other classes of Hamilto- 
nian problems such as the Korteweg-de Vries equations and its generalizations, or nonlinear Klein-Gordon 
equations. The latter also offer possibilities (e.g. in the realm of stability of traveling waves etc.) to consider 
quadratic pencils instead of linear ones, whereby extensions of the Krein matrix method would be partic- 
ularly desirable to develop from a rigorous mathematical viewpoint. These themes are currently under 
study and will be reported in future publications. 
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